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Appraising the projected 10 -18 fractional uncertainty in the optical frequency standards using 
singly ionized ions, we estimate the black-body radiation (BBR) shifts due to the magnetic dipole 
(Ml) and electric quadrupole (E2) multipoles of the magnetic and electric fields, respectively. Mul- 
tipolar scalar polarizabilities are determined for the singly ionized calcium (Ca + ) and strontium 
(Sr + ) ions using the relativistic coupled-cluster method; though the theory can be exercised for 
any single ion clock proposal. The expected energy shifts for the respective clock transitions are 
estimated to be 4.38(3) x 10~ 4 Hz for Ca+ and 9.50(7) x 10~ 5 Hz for Sr+. These shifts are large 
enough and may be prerequisite for the frequency standards to achieve the foreseen 10~ 18 precision 
goal. 

PACS numbers: 31.25.-v, 32.10.Dk, 06.30.Ft, 32.80.-t 



I. INTRODUCTION 

Optical transitions with ultranarrow frequencies in the 
single positively charged ions which use advanced laser 
cooling and trapping techniques are of current interest for 
frequency standards 0-0. Some of the successful single 
ion optical atomic clocks are Hg+ 0,0, Ca+ i> Sr+ 
0,1], A1+ 0,|1], Yb+ QjJ etc., among these the fractional 
uncertainty in Hg + and Al + have already reached 1CU 17 
0| . A new range of experiments are also proposed for 
other ions like Ba+ [HEl, Ra+ [3, In+ [3 El, Yb+ 
16] etc. More precise frequency standards will open-up 
possibilities to study the underlying physics related to 
fundamental constants, probing new elementary physics, 
more importantly it can help in improving the present 
day global positioning svstems and also tracking of deep- 
space probes [H, [l7Hl9| . 

One of the major fortifications to attain smaller frac- 
tional uncertainties in the optical frequency standard 
measurements using ions is the accurate measurement 
of black-body radiation (BBR) shifts. The considered 
standard frequency is shifted from the atomic resonance 
value due to the interaction of the ion with the external 
stray electromagnetic fields present in and around the ex- 
perimental apparatus [2(J. The BBR shift is caused by 
the isotropic field radiated due to finite temperature of 
the apparatus [2014221 ]. 

The dominant contribution to the BBR-induccd energy 
shift is from the electric dipole (El) component of the ra- 
diation field; which have been gauged by many groups for 
a number of ions using the relativistic theories [HI, HH- 
[28j . However, there is absolutely no rigorous estimate 
of the BBR shift due to higher multipoles for any pro- 
posed scheme. Following derivations for El BBR shift in 
|20l - l22| , Porsev and Derevianko have given a generalized 
derivation [25[ to deduce the BBR shifts spawned by any 
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multipole component of the electromagnetic field. 

Two of the proficiency gadgets for the optical frequency 
standards are with a single calcium ion ( 43 Ca + ) |29j 
trapped in a Paul trap and with a strontium ion ( 88 Sr + ) 
confined in an endcap trap 0, H3| . The considered clock 
transitions in these ions are the Si/ 2 —> d 5 / 2 transitions 
operating in the optical regime, the principles are similar 
to the proposed Ba+ [HGj], Ra+ S3 and Yh+ QUE! 
based frequency standards. A major advantage of us- 
ing 43 Ca + ion is that the radiation required for cooling, 
repumping, and clock transition is easily produced by 
non-bulky solid state or diode laser [29]. The reported 
frequency measurements of the transition for frequency 
standard in 88 Sr + have achieved a spectral resolution of 
better than 1.5 Hz I30l - l32l ]. As has been proclaimed, 
these experiments have the dexterity to diminish the rel- 
ative systematic uncertainties to a level of 1CU 17 or below 
7, 13]. In such a scenario, it is compelling to estimate 
the BBR shifts caused by the higher multipoles, espe- 
cially through the Ml and E2 channels for the experi- 
ments comprising the Sx/2 — > d 5 / 2 transitions as in the 
above ions. Such an attempt was made by Porsev and 
Derevianko [25[ for divalent atoms like Mg, Ca, Sr, and 
Yb. In this paper, we extend the work to preview the 
BBR shifts commenced by the Ml and E2 multipoles in 
the ions; particularly for the 43 Ca + and 88 Sr + ions where 
the efficacious experiments are in progress to reduce the 
uncertainties over their previous measurements 0, [29l - 

m. 

The paper is organized as follows: First, we discuss 
the inception of the BBR shifts which embodies the Ml 
and E2 contributions of the radiation field. Then, we 
discuss the method of calculation in the 43 Ca + and 88 Sr + 
ions subsequently presenting the results for the above 
multipolar contributions to the BBR shifts in these ions 
before recapitulating the work. Unless stated otherwise, 
we use the conventional system of atomic units, a.u., in 
which e, m , 47reo and the reduced Planck constant h have 
the numerical value one. 
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II. THEORY 

The interaction Hamiltonian between an electron in an 
atomic system with the external propagating electromag- 
netic field in the Coulomb gauge coupling is given by 



V(r,uj) = — ca • A(r, uj) 

= — c(a ■ e) exp (ik • r), 



(1) 



where a is the Dirac matrix in operator form, uj is the 
angular frequency of the field and k = kk and e are its 
wave vector and polarization direction, respectively. In 
the presence of this interaction, the energy shift that can 
occur for an atomic energy level \^ n ) with energy E n — 
uj n can be approximated to [2lL [22| 



(w n - u m ) 2 - uj 2 



(2) 



Using the multipolar expansion and in terms of the tra- 
ditional multipole moments Q EM , we have [33j 

(a.e)exp(*k.r) = -^^^[Y X LM (k) ■ i] 



LM 



47r(2L + l)(£ + l) ^ 



(3) 

where A = 1 and A = represents the electric and mag- 
netic multipoles, respectively. 

Using the relations B = V x A and E = ^A and 
combining all the above expressions, the energy shift for 
an isotropic field after averaging over u) for all the polar- 
ization and propagation directions is given as [22I [25[ 



SE^ = - 



(aT) 



2L + 1 



2 >Jn 



> 



x F L 



T 



(4) 



where J„ is the angular momentum of state \^ n ) and 



Fl{v) 



L + l 



1 



tt L(2L+ 1)!!(2L- 1)!! 



y 



y-x 



r 2L + l 



1 



-dx. 



(5) 



Here argument y = uj mn /T = (oj n — ui m )/T. The func- 
tion Fl(jj) is a universal function applicable to all the 
atoms with argument y depending on the range of the 
atomic parameters. These functions were first introduced 
by Farley and Wing [22] in the El case and are extended 
in a general form in this work to emanate simpler forms 
for the energy shift expressions. In the limit \y\ » 1 
which corresponds to the case when the transition en- 
ergy is much larger than the temperature (T) is of our 
current interest. 



Substituting values from Eq. (|A4l) in Eq.flU, the BBR 
shift for L = 1 can be expressed as 



1 /8ir 3 a 3 {k B TY 



45(2J„ 
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(6) 

Now at the room temperature (300 K), the BBR shift 
for the El channel can be given as 



8E, 



El 



1 4ir 3 a 3 
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(k B Tfa 
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(831.9V/m) 
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T(K) 
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(7) 



whereas for the Ml channel can be reduced to 
SEr = ~^-(k B T)^ 



2 15 
1 

2 

1 
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(2.77 x 10~ 6 Tesla) 2 



T(K) 
300 



a 



Ml 



(8) 



where a^ 1 , ce^f 1 , (E E1 (uj)) and (-Bf n (w)) are the scalar 
El polarizability, scalar Ml polarizability, the averaged 
El induced electric field and the averaged Ml induced 
magnetic fields, respectively. 

Similarly, substituting values from Eq. (|A5I) in Eq. 
(|4]), the BBR shifts for L = 2 comes out as 

SE (X,2) = _1 ( l6(anf(k B T) e \ ^ | < ^\\Q^\\%n > 
2\ 945(2 J n + 1) J 2^ u mn 

(9) 

which corresponds to the E2 and M2 channels. 

Therefore, the BBR shift at the room temperature for 
the E2 channel is given by 



5E 



E2 



1 8(an) 5 (k B T) & E2 
"2 189(2J„ + 1) 



--(E% 2 (u>))a% 



E2 



1 



(7.2 x 10~ 3 V/m) 5 



T(K) 



1 6 



300 



<\ (10) 



where a^ 2 is the scalar E2 polarizability and (E E2 (lu)) 
is the averaged E2 induced electric field. 

Accumulating all the above expressions, the BBR shifts 
due to both the Ml and E2 channels for an atomic transi- 
tion |\& f) — ¥ l^j) at room temperature can now be given 
in the forms 



--(2.77 x 10~ 6 Tesla) 2 



T(K) 



300 



A/1 



a 



(11) 
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III. METHOD OF CALCULATION 

In order to determine Ml and E2 polarizabilities in 
the considered systems, which have closed-shell configu- 
rations with one valence electron each, we first calculate 
the Dirac-Fock (DF) wave function (|$o)) f° r the cor- 
responding closed-shell configurations, then append the 
valence orbital (n) to define a new reference state (i.e. 
\$n) = a nl^o}) in the Fock space. To obtain the ex- 
act atomic wave functions (ASFs), the core, core- valence 
and valence correlation effects are accounted by defining 
wave operators i7 c , VL cn and Q n , respectively, i.e. the 
exact ASF |^ n ) with a valence orbital n is given by 

= a+O c |$ ) + O cn |$„) + O n |$ n ). (13) 

We adopt the relativistic coupled-cluster (RCC) method 
in the Fock-space to determine ASFs. With T and S n 
representing the core electrons and core electrons with 
the valence electron excitation operators, respectively, 
the above ASF in the RCC framework can be encap- 
sulated in a form [H, [LI M, HI 



(14) 



|*„) = e T {l + S n }|$„>, 



where the ASF for the closed-shell configuration is given 
by Ci c defining |tf„) = ft c |$ > = e T |$ ), ^c„ = e T (™)|$ ) 
and £2 n — e S n \Q n ). We have mentioned T(n) for fi cn 
to emphasize on the fact that the T operator in this case 
excites at least one of the core electrons to the valence 
orbital. All these correlation effects are coupled in the 
course of the ASF determination. 

The equations determining the coupled-cluster ampli- 
tudes and energies are accustomed in compact forms as 



(<S>* \{He T }\<S> a ) = <5o,*A£ c 



(15) 



and 



{K\{HeT}{l + S n }\<S> n ) = 



«|l + 5 n |$ n ) 

($ n \{He T }{l + S n }\$ n ) 
($;|5„,* + S„|$„)A£;„, (16) 



where the notation (* = 1,2) represents for the excited 

hole-particle states, He T denotes the connected terms of 
the Dirac-Coulomb (DC) Hamiltonian (H) with the T 
operators, AE corr and AE V are the correlation energy 
and attachment energy (also equivalent to negative of 
the ionization potential (IP)) of the electron of orbital v, 



respectively. It can be noticed that the reference states 
$o) in Eq. (JTSJ) and |<£„) in Eq. (JTBJl contain different 
number of particles, hence the Hamiltonian used in the 
respective equations describe different number of parti- 
cles in our Fock space representation. We have considered 
contributions only from the singly and doubly excited 
states along with some important valence triple excita- 
tions which are included perturbatively; the approach is 
known as CCSD(T) method in the literature (e.g. see 

To calculate the scalar polarizabilities, it is precedence 
to adopt an approach similar to [1^, [27| or it would be 
prudent to follow a procedure given in (28[. Such ap- 
proaches may be required for achieving better accuracies. 
Accumulating from the previous expressions, the general 
definition of the scalar polarizability is given by 



1 



1 



a 2(A-l) (2£ + l)2J„ + l 



X 



En Ern. 



(17) 



Procuring the derivations in Eq. 
write 



we can now 



a% x = o£* (c) + o£ x (nc) + a£* (n), (18) 

where aS x {c), an x (nc) and a„ x (n) are referred to as 
core, core- valence and valence correlation contributions, 
respectively. 

In the sum-over-states approach as given in Eq. (fT7| . 
it is convenient to determine the low-lying singly excited 
states |\P m ) with respect to the \*f? n ) states that are of our 
interest for both the considered ions. Contributions from 
these states corresponds to the above mentioned valence 
correlation contributions, which are the dominant ones 
compared to the contributions to the part arising from 
higher level excited states that can be estimated from a 
lower order perturbation theory. Contrasting to this con- 
tribution, it is not possible to estimate the core and core- 
valence correlation contributions in a similar procedure. 
As it was stated above, methods describing in [26T - [28j 
would be befitting to account them rigorously. Notwith- 
standing this fact, it will not alter the BBR shifts what 
is the primary intent of the work. This is because of two 
reasons: First, the core correlation effect may be notable 
but in the BBR shift estimation this contribution cancels 
out between the transition states. However, we calculate 
them using the third order relativistic many-body per- 
turbation theory (MBPT(3)) to present the total polar- 
izability results. Secondly, it is observed in the earlier 
works that the core- valence contributions are minuscule 
in the dipole polarizability determinations [H, [24|, [2(| [2?J 
which is further verified below in the present work. So 
it is not necessary to employ a powerful method at the 
cost of heavy computation to determine these tiny con- 
tributions for which we use again the MBPT(3) method 
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to estimate these contributions within the required pre- 
cision level. In this approach, we rewrite Eq. (jTTJ) as 



1 



1 



a 2(X-i) (2L + 1) 2 J„ + 1 



where > is the first order perturbed wave function 

for l 1 ^) due to Qa which is obtained by solving the fol- 
lowing equation 

(H-E n )\¥ r p> = (20) 

where En is the energy correction due to Q\. Some 
of the important diagrams representing this equation for 
the core and core-valcncc correlation effects evaluation 
are shown in Fig. [TJ 



Core-valence correlation evaluating MBPT diagrams 




Core correlation evaluating MBPT diagrams 




FIG. 1: Few important many-body perturbation theory 
(MBPT) diagrams used for the core-valence and core correla- 
tions estimation. 



and 

(k/ || e2 || «i) = < K/ ||C (2) ||^> 

drr^Pfi^P^+Qf^Q^r)), 

(23) 

where P(r) and Q(r) represent the radial parts of the 
large and small components of the single particle Dirac 
orbitals, respectively. The reduced Racah coefficients are 
given by 

M|C (fc) ll«i> = (-l)^ +1/2 V /(2j / + l)(2 J . 1 + l) 

//2 — 1/2) 7r ^ K /' k> (24) 



with 



Tr(l,m,l') = 



1 for I + m + I' 



otherwise. 



(25) 



We have used Gaussian type orbitals (GTOs) to con- 
struct the single particle orbitals for the Dirac-Fock 
(|$o)) wave function calculation. The large and small 
components of the Dirac orbitals in this case are ex- 
pressed as 



and 



d_ 

dr 



(26) 



(27) 



The reduced transition matrix element of a physical 
operator Q\ between |\1//} and |\l/,) in our approach is 
calculated using the expression 



\Qx\ 



(*/ll{i + sJ}Qj{i + ^}ll*i) 



(21) 



where Q\ 



e T 'Q L x e T and M n = <$„|e T 'e T + 
S^e^ e T S n \Q n ) involve two non-truncating series in the 
above expression. We have mentioned concisely the pro- 
cedure for calculating these expressions in the appendix 
in [34||. We evaluate them keeping terms up to fourth 
order in Coulomb interaction. 

The single particle orbital reduced matrix elements for 
the corresponding Ml and E2 cases are taken as 

Ml ml II*) = ^l±M ( _ K/ ||c«||^) 



drriPf^Q^+Qf^P^r)), (22) 



where the summation over k is for the total number of 
GTOs used in each symmetry, and are the nor- 
malization constants for the large and small components, 
respectively, and we use (#• + -) operator to expand 
the small component Dirac orbitals to maintain the ki- 
netic balance condition with its large component. In 
the present calculations, we have considered 9 relativistic 
symmetries (up to g symmetry) and 28 GTOs for each 
symmetry to generate the orbitals. In order to optimise 
the exponents to describe orbitals from various symme- 
tries in a smooth manner, we use the even tempering 
condition 



ah = ao/3 



fc-i 



(28) 



where ao and (3 are two arbitrary parameters that are 
chosen suitably for different symmetries. 

We have considered ao = 7.5 x 10~ 4 for all the sym- 
metries and f3s are optimised to be 2.56, 2.58, 2.61, 2.75 
and 2.83 for s, p, d, f and g orbitals, respectively. For 
the RCC calculations, we have considered excitations up 
to first 16s, 16p, I6d, 14/ and 13<? orbitals while contri- 
butions from other virtual orbitals having large energies 
are neglected. 
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TABLE I: Contributions to the 4si/2 and 3d 5 /2 scalar (a M1 ) 
static polarizabilities in 43 Ca + and their uncertainties in units 
of a 3 ,. The values of the corresponding Ml matrix elements 
for aQo are given in eao- 

Transition Amplitude a An 



TABLE II: Contributions to the 5si/2 and 4d 5 / 2 scalar (a M1 ) 
static polarizabilities in 88 Sr + and their uncertainties in units 
of ag. The values of the corresponding Ml matrix elements 
for aQo are given in eao- 

Transition Amplitude a 



Ot n OL n 



4«i /2 — > 


5si io 

1/ z 


0.0018(3) 


0.4(1) x 10~ 5 


5Si in — > 


6Sl /2 


1.35(2) x 10" 3 


2.79(8) x 10" 6 




6Si/2 


0.0013(2) 


1.7(5) x 10~ 6 




7Sl/2 


1.24(1) x 10~ 3 


1.73(3) x 10~ 6 




3^3/2 


0.0008(1) 


3.1(8) x 10~ 6 




8S1/2 


1.00(1) x 10 -3 


1.0(2) x 10~ 6 




4^3/2 


0.0003(1) 


0.9(7) x 10~ 7 




9Si/2 


1.00(1) x 10~ 3 


9.4(2) x 10~ 7 




5^3/2 


0.0001(1) 


0.8(20) x 10~~ 8 




4^3/2 


5.0(3) x 10" 5 


1.3(2) x 10~ 8 




6^3/2 


0.0003(1) 


0.6(5) x 10~ 7 






1.0(1) x 10 -5 


1.4(3) x 10~ 10 


O-c 






5.0(2) x 10~ 5 




6^3/2 


1.0(1) x 10" 6 


1.1(2) x 10~ 12 








-1.4(3) x 10" 8 




7d 3 / 2 


3.0(2) x 10" 6 


8.8(6) x 10~ 12 








1.0(2) x 10 -8 


Oi c 






1.0(1) x 10~ 3 


C^total 






5.9(3) x 10 -5 








-2(1) x 10~ 8 










Oxtail 






1.0(1) x 10~ 8 


Otn 








Qtotal 






1.0(1) x 10~ 3 


3d 5/2 -4 


3^3/2 
4^3/2 


1.543(5) 
0.004(1) 


-957(6) 
8(4) x 10" 6 


OLn 










5^3/2 


0.004(1) 


8(4) x 10~ 6 


4d 5 / 2 -> 


4d 3/2 


1.545(6) 


-208(2) 




6^3/2 


0.003(1) 


2(2) x 10" 6 




5^3/2 


0.009(2) 


5(2) x 10~ s 




4d 5 /2 


0.004(1) 


8(4) x 10~ 6 






0.006(1) 


1.5(5) x 10~ 5 




5^5/2 


0.009(2) 


3(1) x 10" 5 




7^3/2 


0.004(1) 


7(3) x 10~ 6 




6^5/2 


0.002(1) 
3.2 x 10~ 7 


1(1) x 10~ 6 




5d 5/2 


0.003(1) 


6(4) x 10~ 6 




5ff7/2 


~ 




6d 5 / 2 


0.003(1) 


4(3) x 10~ 6 




6ff7/2 


4.6 x 10" 7 


~ 




7d 5 /2 


0.002(1) 


2(2) x 10~ 6 


Otc 






5.0(2) x 10~ 5 




5ff7/2 


3.2(5) x 10~ 7 


~ 


Qcn 






0.0 




6ff7/2 


3.0(4) x 10~ 7 


~ 


Otail 






-5.0(6) x 10~ 6 


a c 






1.0(1) x 10" 3 


Qtotal 






-957(6) 


OLcn 






-2.5(2) x 10 -7 



Qtail -1-5(1) X 10~ 7 

Ototal -208(2) 



IV. RESULTS AND DISCUSSIONS 



As has been emphasized before, precise estimation of 
the BBR shifts due to Ml and E2 transitions in Ca + 
and Sr + are the focus of this work. The uncertainties 
in these estimations are of two folds: First, the errors 
associated with the considered excitation energies (EEs). 
Second, inaccuracies from the estimated transition ma- 
trix elements. To abate the uncertainties in the evalua- 
tion of the BBR shift for the atomic clock application, 
we use the experimental EEs from the NIST database 
35] for the important singly excited states. Matrix el- 
ements for these states are mentioned explicitly along 
with their uncertainties and respective contributions to 
the the scalar polarizabilities in Tablesfl] UH IIII1 and IIVI 
Contributions from the doubly excited states constructed 
within the configuration space spanned by the considered 
orbitals for our RCC calculations and from the higher 
singly excited states which are not mentioned specifically 
are accounted as at a u using the MBPT(3) method. As 
shown in Tables HI HH IIII1 and IIVI these contributions 
are rather diminutive to be concerned about. The un- 
certainties in these results and for the transition matrix 
elements are estimated taking into account the neglected 



contributions and numerical calculations. 

We first present the static scalar a M1 polarizabilities 
in the Ca + and Sr + ions in Table Q] and Table HH re- 
spectively, for both the ground and d 5 / 2 states. As seen 
in these tables, a M1 s are very small for the correspond- 
ing ground states in both the ions, however they are re- 
ported in this paper for the completeness of the results. 
Contributions from the core correlations are the largest 
in determining these results and they will be cancelled 
while estimating the BBR shifts due to Ml multipole. 
Therefore, the ground state Ml polarizability contribu- 
tions in the considered ions can be neglected. Essentially, 
the contributions to the a M1 polarizabilities in the d 5 / 2 
states are overwhelmingly dominant owing to the very 
small energy gap of their fine structure partner states 
and contributions from all other states are small. The 
Ml matrix elements between the d^/ 2 — d§/ 2 transitions 
in the considered ions were also reported earlier by us 
[3l| and they seem to be very consistent. From our cal- 
culations, we obtain a M1 for the 3^/2 and 4^/2 states 
as -957(6) a.u. and -208(2) a.u. in the Ca+ and Sr+ 
ions, respectively. The uncertainty in the final result is 
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TABLE III: Contributions to the 4s 1( / 2 and 3d 5 / 2 scalar (a E2 ) 
static polarizabilities in 43 Ca + and their uncertainties in units 
of Oo- The values of the corresponding E2 matrix elements 
for Q\ are given in ea\. 

Transition Amplitude a E2 



V O .J 

-> 3d 3/2 


0.12(0] 


212(6) 


4d 3 / 2 


12.51(8) 


121(2) 


5^3/2 


3.89(4) 


9.1(2) 


6^3/2 


5.44(6) 


16.2(4) 


3^5/2 


9.97(6) 


318(3) 


4d 5 /2 


15.30(9) 


181(2) 


5d 5/2 


4.75(5) 


13.6(3) 


6d 5/2 


6.67(6) 


24.2(4) 


a c 




6.15(8) 






0.0 


Otail 




5.36(5) 


Qtotal 




906(5) 


)ther works 




712.91", 871 6 



► 


4S1/2 


9.97(6) 


-106(1) 




5S1/2 


4.99(5) 


9.4(2) 




6S1/2 


1.22(2) 


0.38(1) 




3d 3 / 2 


3.90(4) 


-3657(75) 




4d 3 /2 


4.32(5) 


6.3(1) 




5d 3/2 


1.38(2) 


0.47(1) 




6d 3/2 


1.08(1) 


0.25(5) 




4d 5/2 


8.63(5) 


25.2(3) 




5d 5 / 2 


2.75(3) 


1.87(4) 




6d 5/2 


2.16(2) 


1.02(2) 




5ff7/2 


1.58(2) 


0.56(1) 




6fl7/2 


1.04(1) 


0.210(4) 




5fl9/2 


5.57(5) 


7.0(1) 




6ff9/2 


3.67(4) 


2.62(6) 








6.15(8) 








0.18(2) 


Otail 






-4.29(7) 


Ototal 






-3706(75) 



References: a [3JJ 
b [38] 



obtained by adding the uncertainty from each contribu- 
tion in quadrature. As shown in Table |TT] the uncertainty 
in the d§ / 2 state polarizability value is dominated by the 
uncertainty in the d 3 / 2 — d 5 / 2 contribution. 

We now turn to the a E2 results. Our calculated re- 
sults for the ground and 3d 5 / 2 states in Ca + are given in 
Table Mil Both the results are comparatively large with 
opposite signs. The largest contribution to the ground 
state polarizability comes from the 3d states followed by 
the Ad states. The core contribution is comparatively 
small and the tail contribution to the ground state po- 
larizability is zero due to the absence of occupied d states 
in Ca + . There are two more calculations available in the 
literature on the same using the RCC methods with dif- 
ferent level of approximations [37l l38j . Result reported in 



37] is ab initio and are obtained from a linear response 
theory based calculation. Linear approximation in the 
RCC method is being considered to evaluate the corre- 
sponding transition matrix elements for the estimation of 
the ground state polarizability in [38| using a sum-over- 
states approach like the present work. All the results 
are of same order in magnitude. The ab initio result 
seems to be little lower than the results obtained from 
the sum-over-states approach due to the additional un- 
certainties from the calculated energies that justifies for 
considering the experimental energies in these calcula- 
tions. However, it can be noticed that the differences in 
these results will not alter the BBR shift results which 
can be apparent from the following finding on the E2 
polarizability contributions. To our knowledge, no other 
quadrupole polarizability result is available for the 3d 5 / 2 
state in Ca + to compare with ours. Again, contribution 
from its fine structure partner is the decisive factor for 
the final result followed by a significant contribution from 
the ground state. In Table ITVl we present a E2 results for 
the ground and 4g? 5 / 2 states in Sr + . The magnitudes of 
the ground state result in this ion is larger than Ca + , 
while for the corresponding d 5 / 2 it is other way around. 
There are no results available to the best of our knowl- 
edge to compare our results with them. The trend of the 
contributions from different transitions is almost similar 
for corresponding states in both the ions. We add the un- 
certainty from each contribution in quadrature to obtain 
the final uncertainty in E2 polarizability values. 

Using the above values of the polarizabilities, we ob- 
tain the BBR shift due to the Ml multipole for the 
4s 2 Si/2 — > 3c? 2 £>5/ 2 transition in Ca + at the room 
temperature (300 K) to be 4.38(3) x 10~ 4 Hz. Simi- 
larly, this result comes out to be 9.50(7) x 10~ 5 Hz for 
the 5s 2 Si/2 — > Ad 2 D 5 / 2 transition in Sr + . Contribu- 
tions from the E2 multipole are very small and below the 
uncertainties of the above results and hence can be ne- 
glected for the present purpose of the work. However, the 
reported quadrupole polarizabilities for all the considered 
states may be useful elsewhere. It has to be noticed that 
the BBR shifts due to the El multipole are 0.38(1) Hz 
i|] and 0.22(1) Hz [H for the corresponding transi- 
tions in Ca + and Sr + , respectively. 

V. CONCLUSION 

In summary, we have estimated the black-body ra- 
diation shifts due to the magnetic dipole and electric 
quadrupole multipoles for the 4s 2 S\/ 2 — > 3d 2 D 5 / 2 and 
5s 2 S\/ 2 — > Ad 2 D§/ 2 transitions in the singly ionized cal- 
cium and strontium, respectively. The contribution due 
to the former is the decisive in this case. Nevertheless, the 
reported polarizabilities for the considered states which 
are rarely studied in the above ions may also be useful 
for other purposes. It may be imperative to contemplate 
the reported shifts which are given as 4.38(3) x 10~ 4 Hz 
and 9.50(7) x 10~ 5 Hz in the considered ions to achieve 
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TABLE IV: Contributions to the 5si/ 2 and 4d 5 / 2 scalar (af 2 ) 
static polarizabilities in 88 Sr + and their uncertainties in units 
of Oq. The values of the corresponding E2 matrix elements 
for Q\ are given in ea 2 ,. 



Transition 



Amplitude 



5s 



1/2 





4d 3 /2 


11.25(7) 


382(5) 




5d 3 / 2 


12.87(8) 


136(2) 




6^3/2 


5.00(5) 


16.2(3) 




7^3/2 


3.11(4) 


5.7(1) 




4d 5 /2 


13.91(8) 


572(6) 




5^5/2 


15.64(9) 


201(2) 




6d 5 / 2 


5.97(6) 


23.3(4) 




7d 5 / 2 


3.76(4) 


8.3(1) 


Clc 






14.50(9) 


&cn 






-1.7(2) x 10 


Qtail 






6.35(8) 


Qtotal 






1366(9) 


a„ 



Ad 



5/2 



> 


5S1/2 


13.91(8) 


-191(2) 




6S1/2 


1.31(1) 


0.77(1) 




7S1/2 


1.90(2) 


1.05(2) 




8S1/2 


0.96(1) 


0.231(5) 




9Si/ 2 


0.72(1) 


0.119(4) 




4d 3 / 2 


6.08(6) 


-1930(38) 




5^3/2 


5.63(6) 


12.1(3) 




6d 3 / 2 


1.83(2) 


0.92(2) 




7c?3/2 


1.14(1) 


0.315(6) 




5c?5/2 


11.17(7) 


47.4(6) 




6c?5/2 


3.65(5) 


3.7(1) 




7d 5 / 2 


2.28(3) 


1.27(3) 




5^7/ 2 


2.73(3) 


1.76(4) 




657/2 


2.21(2) 


1.11(2) 




5fl9/2 


9.64(7) 


22.0(3) 




659/2 


7.80(6) 


13.7(2) 








14.50(9) 








0.24(3) 


Qtail 






-4.83(5) 


Ototal 






-2005(38) 



the 10 18 precision uncertainty in the proposed clock ex- 
periments. 
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Appendix A: Farley & Wing's functions 

With the aid \y\ >> 1, the following the expression 
1 L + l 



F L (y) 



n L(2L + 1)\\(2L - 1)!! 

f°° ( 1 1 \ x 2L+1 

I l^ + ^)^d x ,A1) 



^y + x y — x J e x — 1 
gives for L = 1 as 

2 f°° ( 1 1 



Fi{y) 



3ttJo \y + x y — x J e x — 1 

2 (2 f°° x 3 , 
-ax 



-dx 



3tt V y Jo e x - 1 y A J e x - 1 



-dx . 



(A2) 



Further by using the definite integral value 



where £?2« is the Bernoulli number, Fi(y) reduces to 

4tt 4 



Fi(y) 



45y' 



(A4) 



be 



Similarly for L = 2, the above expression turns out to 



1 1 \ X 5 

H I — -dx 



y + x y — x J e x — I 



1% 7o e x — 1 y 3 Jo e x — I 
8tt 5 



945y 



(A5) 



Appendix B: Square of the matrix element 

In our approach, we write 

|*„) - otn c |$o>+ncn|*»>+fin|*»>- (Bl) 

With this expression, the square of the matrix element 
of any arbitrary operator O can be expressed as 

(<f n \0\* m ) 2 = (* n |0|* m > <* m |0|*„) 
= ($o|^OO ro Ot,OO c |$ ) 

+ ($ \nlon cm ni m on c \$ ) 
+ ($ n \nl n on c nlon cn \$ n ) 
+ ($ n \nion c nlon n \$ n ) 
+ ($ n |fitOfi cm fit m on n |$ n ) 

+ ($ n |fitOfi rn fit,Ofi n |* n ) , (B2) 
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where we have facilitated the generalized Wick's theo- 
rem to derive these terms and assumed all the operators 
are in normal order form and only the connected terms 
are survived. For simplicity the first two terms, the sec- 
ond term and the last three terms are categorized into 
core (c), core- valence (cn) and valence (n) contributions, 



respectively; i.e. in an abbreviate form it is given as 

(*„|o|* ro ) 2 = (y n \o\y m )l + (y n \o\* m ) 2 cn 

+ (*n\0\* m f n . (B3) 
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